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ABSTRACT 


An  experimental  apparatus  was  designed,  constructed  and 
instrumented  in  an  effort  to  systematically  and  carefully 
study  the  condensation  heat-transf er  coefficient  on  a  single, 
horizontal  tube.  A  smooth,  thick-walled  copper  tube  of 
length  133.5  mm,  with  an  outside  diameter  of  15.9  mm  and 
an  inside  diameter  of  12.7  mm  was  instrumented  with  six 
wall  thermocouples.  The  temperature  rise  across  the  test 
section  was  measured  accurately  using  quartz  crystal 
thermometers.  The  inside  heat-transf er  coefficient  was 
determined  using  the  Sieder-Tate  correlation  with  leading 
coefficient  of  0.029.  Initial  steam  side  data  were  taken 
at  atmospheric  pressure  to  test  the  data  acquisition/ 
reduction  computer  programs. 
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I.  INTRODUCTION 

A.  SIGNIFICANCE 

Experimental  studies  have  shown  that  enhancement  techniques 
applied  on  the  steam  side  of  horizontal  condenser  tubes  may 
significantly  reduce  the  thermal  resistance  and  thus  reduce 
the  overall  size  and  weight  of  the  condenser.  Several  work¬ 
shops  and  conferences  on  heat  exchangers  have  addressed 
specific  areas  in  need  of  research  {Refs.  1,2,3).  Recently, 
there  have  been  many  investigations  into  possible  alterna¬ 
tives  for  enhanced  heat-transf er  surfaces  in  the  design  of 
heat  exchangers  {Refs.  2, 4, 5, 6, 7).  Webb  {Ref.  8}  in  1980 
published  an  excellent  review  of  enhancement  methods  as 
applied  to  condensers. 

Enhancement  techniques  can  generally  be  divided  into 
two  major  categories:  water-side  enhancement  and  steam-side 
enhancement.  Water-side  enhancement  is  covered  extensively 
by  Bergles  &  Jensen  {Ref.  9).  Various  steam-side  enhance¬ 
ments  include  low  integral  fins,  surface  coatings,  fluted 
tubes  and  roped  tubes. 

The  goal  of  any  enhancement  technique  is  to  lower  the 
thermal  resistance  to  the  point  where  the  technique  becomes 
economically  feasible.  Althouth  Bergles  &  Jensen  {Ref.  9} 
show  that  extensive  reduction  in  thermal  resistance  can  be 
obtained  on  the  water-side,  it  is  felt  that  a  significant 
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portion  of  the  overall  thermal  resistance  can  be  reduced  by 
enhancement  techniques  on  the  steam-side.  Both  methods  must 
be  examined  independently  because  the  best  enhancement  for 
one  side  may  not  be  the  best  for  both  sides. 

Several  research  efforts  have  been  undertaken  at  the 
Naval  Postgraduate  School  to  model  steam  surface  condensers 
and  analyze  their  thermal  performance  numerically  {Refs. 
10,11,12}.  Based  on  numerous  computer  runs  used  to  analyze 
different  combinations  of  both  inside  and  outside  enhancements, 
it  is  felt  that  further  experimental  research  should  be 
conductdd  on  steam-side  enhancement  techniques. 

B.  OBJECTIVE 

The  objective  of  this  research  effort  was  to  obtain 
smooth-tube  baseline  data  necessary  to  act  as  the  standard 
for  future  enhancement  techniques.  Both  the  inside  and 
outside  heat-transfer  coefficients  had  to  be  identified 
and  all  experimental  methods  used  to  obtain  and  reduce  the 
recorded  data  had  to  be  well  documented  so  that  future 
testing  using  enhanced  tubes  would  be  performed  under  identical 
conditions. 
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II.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 


A.  OVERVIEW  OF  SYSTEM 

The  experimental  apparatus  used  in  this  investigation,  as 
described  in  Reference  13.  was  constructed  of  stainless  steel 
and  152.4  nm  (6  in)  diameter  Corning  Pyrex  glass  piping.  An 
overall  schematic  of  the  system  is  shown  in  Figure  2.1  and 
a  photograph  is  shown  in  Figure  2.2.  A  boiler  section 
304.8  mm  (12  in)  in  diameter  and  487.7  mm  (19.2  in)  in 
height,  with  a  304.8  mm  (12  in)  to  152.4  mm  (6  in)  reducing 
section,  provides  steam  for  the  system.  Ten  Watlow 
4000-Watt,  480-Volt  stainless  steel  immersion  heaters  are 
used  to  provide  input  power.  Each  heater  is  15.9  mm  (5/8  in) 
diameter  and  279.4  mm  (11  in)  long.  All  heaters  are 
connected  in  parallel. 

Stean  travels  up  the  glass  piping  to  a  height  of  3.1  m 
(10  ft)  and  makes  a  180-degree  bend  followed  by  a  1.52  m 
(5  ft)  straightening  section  prior  to  passing  through  the 
condenser  test  section.  Two  helically-wound,  water  cooled 
coils  of  9.5  mm  (3/8  in)  diameter  copper  tubing  act  as  a 
dump  condenser  to  remove  all  excess  steam.  A  25.4  mm  (1  in) 
diameter  stainless  steel  tube  returns  condensate  by  gravity 


B.  CONDENSER  TEST  SECTION 


The  condenser  test  section  is  made  of  stainless  steel, 

152.4  mm  (6  in)  in  diameter  and  457.2  mm  (18  in)  in  length. 
Figure  2.3  shows  a  schematic  and  Figure  2.4  shows  a  photo¬ 
graph  of  the  test  section.  Steam  is  condensed  on  a  single 
horizontal  tube.  Two  types  of  smooth  test  tubes  were 
constructed  of  high-grade  copper  to  provide  baseline  data  for 
inside  and  outside  heat-transfer  coefficients.  An  instrumented 
tube,  constructed  as  shown  in  Figure  2.5,  and  an  uninstrumented 
tube  of  the  same  dimensions  were  used  to  obtain  initial  data. 
The  uninstrumented  tube  was  constructed  identical  to  the 
instrumented  tube  but  there  were  no  channels  machined  for 
thermocouple  installation.  The  active  condensing  length 
was  133.3  mm  (5.25  in)  with  a  suitable  correction  applied 
for  the  "fin"  effect  created  by  both  ends  of  the  tube  not 
being  in  direct  contact  with  the  steam  on  the  outside, 
but  cooled  by  water  on  the  inside.  The  fin-effect  calculations 
are  described  in  Appendix  A. 

The  test  tube  is  held  in  place  by  nylon  holding  rings  as 
described  in  Reference  13.  An  additional  0-ring  was  placed 
inside  the  inlet  section  to  prevent  cooling  water  from 
leaking  into  the  grooves  of  the  instrumented  test  tube. 

Nylon  was  used  to  minimize  axial  conduction  from  the  test 
section.  A  nylon  mixing  chamber  described  in  Reference  13 
was  placed  in  the  outlet  side  prior  to  measuring  the  temper¬ 
ature  to  ensure  thorough  mixing  of  the  coolant.  A  114.3  mm 


(4.5  in)  diameter  double  glass  window,  heated  by  hot  air 
to  prevent  fogging,  was  used  to  provide  visual  observation 
of  the  condensing  process. 

A  12.7  mm  (0.5  in)  port  at  the  dump  condenser  section 
was  used  as  an  overpressure  relief  line  and  boiler  fill 
line.  An  in-line  0.108  MPa  (1  psig)  relief  valve  was  used 
to  allow  gases  to  escape  while  warming  up  the  system. 

A  12.7  mm  (0.5  in)  vacuum  valve  was  placed  in  a  branch  line 
off  of  the  relief  line  to  allow  filling  of  the  boiler  as 
described  in  Appendix  B.  All  connections  were  assembled 
with  Swagelok  fittings  sealed  with  Teflon  ferrules  to  ensure 
leak  tighness. 

1 •  Cooling  Water  System 

A  closed  loop  cooling  water  system  was  designed  and 

constructed  to  circulate  cooling  water  through  the  test 

section.  A  0.001  cu-m  (2.8  cuft)  stainless-steel  tank  was 

used  as  a  reservoir  which  fed  the  suction  end  of  a  close 

3 

coupled  centrifugal  pump  rated  at  0.0022  m  /s  (11  GPM)  at 
12.19  m  (40  ft)  of  head.  A  19.1  mm  (3/4  in)  diameter  bronze 
needle  valve  was  used  to  regulate  the  flow  which  was 
monitored  by  a  rota  eter  with  a  full-scale  range  of  6.94 
E-4  m  /s  (11GPM).  See  Figure  2.6  for  details.  The  maximum 
flow  rate  through  the  test  tube  was  4*57  m/s  (15  fps)  at 
an  inlet  temperature  of  15.6°C  (60°F).  The  centrifugal 
pump  needle  valve  arrangement  provided  a  flow  rate  accuracy 
of  0.5  percent  on  a  scale  with  100  divisions.  After  leaving 
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the  condenser  test  section,  the  water  wa~  cooled  back  down 
to  the  15.6°c  (60°F)  bath  temperature  by  means  of  a  11.73  kW 
( 40, 000-BTU/Hr )  direct  expansion  heat  exchanger.  The  heat- 
exchanger  refrigerant  dissipates  the  heat  to  the  atmosphere 
through  a  10.55  kW  (3  ton)  air-conditioning  unit.  The 
chiller  section  is  illustrated  in  Figure  2. 7. 

The  closed  loop  cooling  water  system  was  chosen  in 
order  to  more  closely  monitor  the  water  flow  rate  and  the 
inlet  temperature  to  the  test  condenser  tube.  The  reservoir 
is  easily  filled  by  use  of  the  installed  fill  line  and 
drained  by  use  of  a  siphon  tube. 

2.  Test  Condenser  Tubes 

Initial  smooth  tube  data  were  taken  using  an  instru¬ 
mented  copper  tube  as  illustrated  in  Figure  2.8.  Construction 
of  the  instrumented  tube  was  carefully  performed  in  accordance 
with  Figure  2.5.  The  instrumented  tube  must  be  installed 
according  to  the  instructions  as  outlined  in  Appendix  D  to 
prevent  damage  to  thermocouple  wires  and  leakage  from  the 
system.  Each  wall  thermocouple  should  be  checked  for  proper 
operation  prior  to  taking  any  data. 

C.  DUMP  CONDENSER  SYSTEM 

A  dump  condenser  is  installed  to  condense  excess  steam 
from  the  test  section.  Cooling  water  is  supplied  from  tap 
water  through  two  parallel,  helically-wound,  9.5  mm  (3/8  in) 
diameter,  copper  coils.  One  smaller  diameter  coil  is  placed 
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inside  the  other,  while  both  are  457.2  mm  (18  in)  high.  The 
tap  water  flow  rate  is  controlled  by  a  19.05  mm  (3/4  in) 
diameter  bronze  needle  valve  and  is  monitored  using  a 
rotameter.  System  internal  pressure  can  be  controlled  by 
reducing  the  flow  of  cooling  water  to  the  dump  condenser. 

For  a  given  boiler  power  input,  a  dump  condenser  rotameter 
reading  will  correspond  to  a  system  pressure. 

D.  VACUUM  SYSTEM 

Two  operational  vacuum  systems  are  available  to  maintain 
sub-atmospheric  conditions  within  the  apparatus;  an  air 
ejector,  powered  by  1.1  MPa  (160  psig)  air  supplied  by  the 
Mechanical  Engineering  Department  reciprocating  air  corn- 
pressure,  and  a  portable  mechanical  vacuum  pump  (Figure  2.6). 
Either  system  may  be  used  to  initially  draw  down  the  pressure 
but  only  the  air  ejector  may  be  used  while  the  system  in  in 
operation  due  to  the  lack  of  a  nitrogen  cold  trap  on  the 
mechanical  vacuum  pump.  The  vacuum  system  inlet  valve  can 
be  used  to  throttle  the  air  ejector  to  maintain  any  system 
pressure  above  13.8  kPa  (2  psia).  This  system  is  also 
useful  in  removing  noncondensable  gases  from  the  apparatus. 
The  sequential  drawing  down  of  vacuum  on  the  system  while 
steaming,  shutting  the  vacuum  inlet  valve  and  raising  the 
pressure  to  about  13.8  kPa  (2  psia)  by  adjusting  the  dump 
condenser  throttle  valve  was  found  to  be  an  effective 
method  of  removing  noncondensable  gases.  This  procedure 
should  be  repeated  about  three  times. 
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The  system  will  not  maintain  vacuum  due  to  cracked  silver 
soldered  connections  at  the  base  of  the  dump  condenser 
section.  The  copper  tubing  used  as  a  dump  condenser  pene¬ 
trates  the  stainless  steel  base  and  was  sealed  with  silver 
solder  which  cracked,  producing  small  leaks. 

E.  INSTRUMENTATION 

1 .  Heater  Power  Control 

Power  input  to  the  apparatus  is  regulated  by  a 
voltage  sensing  circuit.  Line  voltage  (440-VAC)  is  fed  into 
a  differential  input  precession  voltage  attenuator  which 
divides  the  line  voltage  by  one-hundred.  This  stepped- 
down  voltage  passes  through  a  True  Root  Mean  Square  converter 
stage  on  which  the  integrated  period  has  been  reduced  to 
about  1  ms.  The  output  of  the  TRMS  converter  is  buffed  and 
compared  to  a  reference  voltage  provided  by  a  potentiometer 
mounted  on  the  front  panel.-  See  Figure  2.9  for  a  general 
view  of  the  control  panel.  The  comparator  output  is  fed 
to  the  control  input  of  a  Halmar  silicon-controlled  rectifier 
power  supply  which  provides  the  voltage  applied  to  the  heater 
elements.  The  buffered  output  of  the  TRMS  converter  is 
amplified  and  filtered  to  provide  a  normalized  voltage 
output  proportional  to  the  actual  voltage  output  of  the 
power  supply  which  is  monitored  by  the  data  acquisition/ 
reduction  system.  A  line  diagram  is  provided  in  Figure 
2.10. 
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2 .  Pressure  Measurement 


A  U-Tube  mercury  manometer,  calibrated  in  millimeters, 
is  installed  to  measure  system  pressure.  The  pressure  tap 
is  located  at  the  test  section  and  is  tilted  upward  to  avoid 
steam  from  collecting  in  the  manometer  (Figure  2.4).  A 
vacuum  valve  is  installed  in  the  pressure  line  to  prevent 
condensation  from  accumulating  in  the  manometer.  A  small 
amount  of  water  eventually  collects  in  the  manometer,  but 
this  discrepancy  is  taken  into  account  by  the  data  acquisition/ 
reduction  software  with  manual  input  to  the  keyboard  of 
mercury  height  followed  by  water  height. 

3 •  Temperature  Measurement 

Various  temperatures  are  monitored  throughout  the 
system  to  include:  1.  test  section  cooling  water  inlet 
and  outlet,  2.  steam,  3.  condensate  return,  4.  ambient, 

5.  test  tube  wall  and  6.  cooling  water  reservoir.  The 
temperature  rise  across  the  test  section  is  measured  by  a 
Hewlett-Packard  quarts  crystal  thermometer  accurate  to 
0.04°C.  The  location  of  the  temperature  probes  are  shown  in 
Figure  2.4.  The  cooling  water  reservoir  temperature  is 
monitored  by  a  direct  readout  alcohol-in-glass  thermometer. 

All  other  temperature  measurement  is  accomplished  by 
copper-constantan  thermocouples:  two  for  the  steam,  six  for 
the  test  tube  wall  and  one  for  the  condensate  return.  Each 
temperature  measurement,  except  for  the  cooling  water 
reservoir,  is  read  directly  by  a  Hewlett-Packard  3497A 
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data  acquisition  system,  which  is  controlled  by  a  Hewlett- 
Packard  9826  computer. 

During  operation,  the  six  wall  thermocouples 
experienced  temperature  fluctuations  on  the  order  of  -  10  qV. 
Therefore,  each  thermocouple  was  scanned  for  ten  seconds 
and  the  ten  readings  were  averaged  to  obtain  a  more  accurate 
measurement.  All  other  thermocouples  maintained  a  steady 
readout  with  only  1.0  qV  fluctuation.  The  quartz  crystal 
thermometer  appeared  to  operate  very  well  with  no  fluctuations. 

Two  of  the  thermocouples  were  calibrated.  One 
thermocouple  was  made  from  the  beginning  of  a  spool  of 
copper-constantan  wire  and  one  was  made  from  the  end. 

All  other  thermocouples  were  fabricated  from  wire  on  this 
spool.  It  was  assumed  that  the  properties  of  the  copper 
and  constantan  do  not  change  along  a  given  section  of  wire 
for  any  given  spool.  Both  thermocouples  were  calibrated 
by  the  method  described  in  Appendix  E. 
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Figure  2.1  Schematic  of  Test  Apparatus 


Steam 

Temperature 


Control 


Cooling  Water 
Reservoir 


Figure  2.6  Photograph  of  Cooling 
Water  System 
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Figure  2.9  Photograph  of  Control  Panel 


III.  SYSTEM  OPERATION 


After  a  thorough  cleaning  using  the  procedures  described 
in  Appendix  F,  a  test  tube  may  be  installed  by  the  steps 
outlined  in  Appendix  D.  Each  system  start-up  should  be 
accomplished  by  the  procedures  described  in  Appendix  B. 

Once  the  system  has  reached  steady-state  conditions  (about 
30  minutes),  data  taking  for  a  given  test  tube  may  commence. 
During  each  data  run,  a  frequent  check  of  the  test  tube 
should  be  made  to  ensure  the  desired  condensation  mode: 
either  dropwise  or  filmwise.  When  the  tube  cleaning  procedures 
of  Appendix  F  are  utilized,  filmwise  condensation  conditions 
will  last  for  about  four  hours,  followed  by  mixed  filmwise 
and  dropwise  conditions.  Each  data  set  consisted  of  eight 
cooling  water  velocities  ranging  from  1.37  m/s  to  3.84  m/s 
through  the  test  section.  At  a  given  velocity  setting,  the 
flow  was  allowed  to  steady  out  for  two  minutes  before 
recording  temperatures.  The  computer  signalled  by  both 
audio  and  visual  means  when  the  system  was  ready  for  a  new 
data  point.  An  entire  data  set  was  recorded  in  40  minutes. 
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IV.  DATA  ACQUISITION /REDUCTION 


A.  DATA  ACQUISITION  AND  STORAGE 

A  Hewlett-Packard  (HP)  3054A  automatic  data  acquisition/ 
control  system  was  used  to  monitor  temperatures  from  the 
thermocouples  and  the  quartz  crystal  thermometers.  Figure  4.1 
shows  a  photograph  of  the  system.  This  system  included  an 
HP  3497A  data  acquisition/control  scanner  and  an  HP  3456A 
digital  voltmeter  with  a  resolution  of  100-nanovolts.  An 
HP  9826a  computer,  interfaced  with  the  data  acquisition/ 
control  system,  was  used  as  a  controlling  unit  and  provided 
storage  for  data.  Information  was  entered  through  the  key¬ 
board  to  prompt  the  data  acquisition/control  unit  to  auto¬ 
matically  scan  each  channel.  Channel  assignments  are  listed 
in  Table  1 .  These  raw  data  were  then  transferred  to  a 
computer  disk  under  a  user-specified  file  name  for  later 
reduction  to  usable  form.  The  ability  to  store  raw  data 
directly  enabled  these  data  to  be  reduced  at  any  time  and 
allowed  flexibility  for  changes  to  data-reduction  software. 

B.  DATA  REDUCTION 

Following  data  acquisition  for  each  data  point,  results 
were  computed  according  to  the  stepwise  procedure  outlined 
in  the  next  section,  and  were  then  printed  on  an  HP  6771G 
thermal  printer.  Appendix  G  shows  an  example  of  a  represen¬ 
tative  data  run.  Results  were  also  stored  in  a  user-specified 
plot  data  file  for  subsequent  plotting  using  the  program  PLOT. 
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C.  STEPWISE  SOLUTION  PROCEDURE 

1 .  Program  SIEDER 

a.  Compute  average  cooling  water  temperature. 

b.  Compute  average  wall  temperature. 

c.  Compute  cooling  water  velocity. 

d.  Compute  mass  flow  rate  of  cooling  water. 

e.  Compute  heat  transferred  to  cooling  water. 

f.  Compute  average  inside  wall  temperature. 

g.  Compute  LMTD. 

h.  Assume  fin  efficiency. 

i.  Compute  inside  heat- transfer  coefficient. 

j.  Compute  Nusselt  number. 

k.  Calculate  fin  efficiency  using  previously- 
calculated,  heat-transfer  coefficient. 

l.  Recalculate  heat-transf er  coefficient. 
Recalculate  Nusselt  number  and  compare  to  Nusselt 

number  found  in  j  and  iterate  if  not  within  Figure  5.1 

shows  an  example  of  plotted  data  and  Appendix  G  shows  a 
listing  of  PLOT. 

2.  Program  NSFDRP 

a.  Compute  average  cooling  water  temperature. 

b.  Compute  cooling  water  velocity. 

c.  Compute  mass  flow  rate  of  cooling  water. 

d.  Compute  heat  transferred  to  cooling  water. 

e.  Compute  log  mean  temperature  difference. 

f.  Compute  overall  heat- transfer  coefficient  based 
on  outside  area. 


g.  Compute  wall  resistance  based  on  outside  area. 

h.  Compute  Reynolds  number  of  cooling  water. 

i.  Assume  fin  efficiency. 

j.  Compute  inside  heat- transfer  coefficient. 

k.  Calculate  cooling  water  temperature  rise. 

l.  Calculate  fin  efficiency  using  previously- 
calculated,  heat-transf er  coefficient. 

m.  Recalculate  heat-transf er  coefficient  and  iterate 
if  not  within  'll!)  of  the  value  found  in  h. 

n.  Calculate  condensing  heat- trans f er  coefficient 
from  Uq,  Rw  and  . 


TABLE  I 

HP  3497A  Channel  Assignments 


Channel 

Assignment 

20 

steam 

21 

steam 

22 

condensate 

23 

room 

24 

room 

25 

wall 

26 

wall 

27 

wall 

28 

wall 

29 

wall 

30 

wall 
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Figure  4-1  Photograph  of  Data  Acquisition/Reduction  System 


V.  RESULTS  AND  DISCUSSIONS 

A.  INSIDE  HEAT-TRANSFER  COEFFICIENT 

Experimental  data  were  taken  at  atmospheric  pressure  by 
the  procedures  outlined  in  Chapter  III.  Twelve  data  sets 
were  taken  to  determine  the  inside  heat- transfer  coefficient. 
The  first  four  data  sets  were  utilized  to  test  the  data 
acquisition/reduction  program  SIEDER  and  the  middle  five 
data  sets  experienced  mixed  filwise  and  dropwise  condensa¬ 
tion  conditions.  Data  sets  10-12  showed  filmwise  condensa¬ 
tion  throughout  each  of  the  data  collection  sets.  Figure  5.1 
is  a  plot  of  data  taken  during  run  number  10.  Similar  results 
were  obtained  for  data  runs  11-12.  The  analysis  of  all 
data  showed  a  Sieder-Tate  coefficient  of  0.029-0.001  to  be 
representative  for  the  experimental  apparatus  (Table  II). 

The  Sieder-Tate  correlation  for  fully  developed  turbulent 
flow  has  a  leading  coefficient  of  0.027  (Ref.  1 6 } .  Rose 

{Ref.  14)  shows  a  leading  coefficient  of  0.030  with  an 
L/D  ratio  of  17.  which  is  in  good  agreement  with  the  results 
of  this  thesis.  The  L/D  ratio  for  the  experimental  apparatus 
used  during  this  thesis  was  18,  well  below  that  for  fully 
developed  flow  (Ref.  16};  therefore  the  inside  heat- transfer 
coefficient  was  correspondingly  higher. 
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TABLE  II 


Sieder-Tate  Coefficients 


File 

Ci 

10 

.0291 

1 1 

.0287 

12 

.0287 

B.  OUTSIDE  HEAT-TRANSFER  COEFFICIENT 

Experimental  data  were  taken  at  atmospheric  pressure  by 
the  procedures  outlined  in  Chapter  III.  Four  data  sets  were 
recorded  to  detemine  the  outside  heat- transfer  coefficient 
using  the  program  NSFDRP  (Appendix  G)  to  reduce  the  raw 
data.  Figure  5.2  shows  a  plot  of  the  four  data  sets  along 
with  a  curve  showing  the  corresponding  Nusselt  prediction. 
All  of  the  data  lie  above  the  theoretical  curve.  This 
discrepancy  has  two  possible  explanations.  Either  dropwise 
condensation  and/or  vapor  shear  associated  with  a  vapor 
velocity  of  3  m/s  could  increase  the  condensation  heat 
transfer  coefficient.  The  effect  of  vapor  shear  has  been 
studied  by  many  investigators  including  Rose  as  described 
in  a  very  recent  report  {Ref.  18}.  Rose  reported  that  when 
vapor  velocity  was  increased  to  7  m/s  for  refrigerant-113, 
the  condensing  coefficient  increased  by  20$.  The  presence 
of  dropwise  condensation  significantly  raises  the  heat- 
transfer  coefficient,  sometimes  as  much  as  ten  times 
{Ref.  17}. 
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A  large  amount  of  scatter  occurred  in  the  recorded  data. 
This  situation  could  have  occurred  due  to  noncondensable 
gases  and  the  inability  to  regulate  system  pressure.  During 
normal  operation,  the  Mechanical  Engineering  Department 
air  compressor  (160  psi)  was  used  to  operate  the  air 
ejector.  During  the  condensation  heat  transfer  coefficient 
runs,  this  air  compressor  was  inoperative;  therefore,  house 
air  (50  psi)  was  used.  This  source  of  air  had  a  very  low 
capacity  and  fluctuated  greatly  during  each  data  run, 
perhaps  causing  varying  amounts  of  air  to  be  in  the  vicinity 
of  the  test  tube. 

Due  to  limited  and  lack  of  reliable  data,  very  little 
can  be  concluded  for  the  results  thus  far  for  the  condensing 
coefficient.  More  data  is  necessary  to  reach  any  final 
conclusions . 
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Figure  5*1  Sieder-Tate  Paramater  Plot 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

1.  An  experimental  apparatus  that  allows  the  study  of 
condensation  heat  transfer  of  steam  on  a  single,  horizontal 
tube  has  been  designed,  constructed  and  instrumented. 

2.  The  water-side  heat- transfer  coefficient  was 
established  using  the  Sieder-T&te  correlation  with  a  leading 
coefficient  of  0.029  -  0.001. 

3.  An  analysis  scheme  to  establish  the  condensing  heat- 
transfer  coefficient  has  been  outlined,  and  data  acquisition/ 
reduction  programs  have  been  written  and  tested. 

B.  RECOMMENDATIONS 

1.  To  prevent  leaking,  install  Swagelok  fittings  with 
Teflon  ferrules  into  the  dump  condenser  lines  penetrating 
the  stainless  base  of  the  dump  condenser  section. 

2.  Rig  suspension  cables  to  remove  the  dump  condenser 
section  by  attaching  them  in  the  ceiling  using  metal  studs. 

3.  Install  a  thermopile  at  the  outlet  of  the  test 

condenser  tube  as  a  second  means  of  measuring  the  cooling 

water  temperature  rise  across  the  test  section. 

/ 

4.  Use  a  commutating  bridge  to  accurately  determine 
boiler  heater  resistance  over  the  full  range  of  power 
settings . 
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5.  Install  a  nitrogen  cold  trap  on  the  mechanical  vacuum 
pump. 

6.  Install  a  permanent  drain  line  on  the  cooling  water 
reservoir  to  replace  the  temporary  siphon  tube. 

7.  Replace  the  extension  wire  from  the  ice  bath  to  the 
data  acquisition  system  with  thermocouple-quality  copper 
wire . 

8.  After  recommendation  No.  1  above  is  completed,  check 
the  system  for  leak  tightness  under  vacuum  conditions. 

9-  Complete  experimentation  to  determine  the  condensing 
heat-transfer  coefficient  at  1  atmosphere  as  well  as  under 
vacuum  conditions. 
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APPENDIX  A 


SAMPLE  CALCULATIONS 

1 .  INSIDE  HEAT-TRANSFER  COEFFICIENT 

A  set  of  calculations  was  performed  to  verify  data  reduc¬ 
tion  program  SIEDER.  One  data  point  was  selected  arbitrarily 
from  data  set  10. 

a.  Given  the  following  inlet  and  outlet  temperatures, 
compute  average  cooling  water  temperature. 

T1  =  23 . 94°C 
T2  =  26.32°C 

Tb  =  (23.94  +  26.32)/2  =  25.13°C 

b.  Compute  average  wall  temperature. 

T  =  (  T  .  +T  -  +  T  _  +  T  ,+T  c  +  T  /■)  /6  =  62.64°C 
w  v  wl  w2  w3  w4  w5  wo' 

c.  Compute  water  velocity.- 

Flow  meter  reading  =  25$  and  from  the  calibration 

curve,  V  =  4.5  ft/s  or  V  =  1.37  m/s. 
w  w 

From  Reference  16,  Table  A. 6  at  300  K  the  following 
properties  are  found: 

Pr  =  5.83 

Uj  =  855E-6  N-s/m2 

C  =  4.179  kJ/Kg-K 
P 

Kf  =  61 3E-3  W/m-K 
vf  =  1.003E-3  m3/kg 

d.  Compute  mass  flow  rate  of  cooling  water. 

A  =  iTDi2/4 

m  =  V  A  If  -  0.173  kg/s 
w 
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e.  Compute  heat  transferred  to  cooling  water. 

Q  -  m  Cp  (T1-T2) 

Q  =  1721  W 

f.  Compute  average  inside  wall  temperature. 

T  =  Q  In  (D  /D.)  /2  i  k.  L  Reference  17 

w  01  f 

T  =  2.1 6°C 
w 

T  .  T  -  0.5  AT 
Wl  =  w  w 

T  .  =  61 . 56°C 

Wl 

g.  Compute  log  mean  temperature  difference. 

LMTD  =  (T2-T1  )/ln  [(Twi-T1  )/(Tw.-T2)] 

h.  A  fin  efficiency  is  assumed  to  be  20$  for  first 
interaction . 

i.  Compute  heat-transfer  coefficient. 
hi  =  Q/t r  D  (L+I^  «n1+L2*n2)  LMTD 

hi  =  7772  W/m2-K 

j.  Compute  Nusselt  number. 

Nu  =  hiDi/kf> 

Nu  =  1 61  -  first  iteration 

k.  Calculate  fin  efficiency. 

P1  =  2tt  DQ  =  0.1197  m 

P„  =  2tt  D  =  0 . 1 00  m 
2  r 

A,  =  it  D  (D  -D.  )  =  0.00038  m2 
1  0  0  1 

A0  =  7 r  D  (D  -D.  )  =  1.6  E-4  m2 
r  r  l 


P^/k^A^ 

m1  =  79.74 

m2  "  ^i  P2//jcfA2 
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m2  =  112.32 

n1  =  tanh  (m1  )/m1  L 

n1  =  .2079 
n2  =  .2913 

l.  Recalculate  heat- transfer  coefficient. 
hj_  =  7589.9  W/m2-K 

m.  Recalculate  Nusselt  number  and  compare  to  previous 
calculated  value. 

Hi  =  157.24 

This  value  for  the  Nusselt  number  agrees  within  0.3$  of 
the  computer  generated  value  (Nuc  =  156.7). 

2.  CONDENSATION  HEAT-TRANSFER  COEFFICIENT 

A  set  of  calculations  was  performed  to  verify  the  data 
reduction  program  NSFDRP.  One  data  point  was  selected  from 
data  run  1  . 

a.  Given  the  following  inlet  and  outlet  cooling  water 
temperatures,  compute  average  water  temperature. 

T1  =  30. 06° C 
T2  =  32.53°C 
Tb  =  (Tn  +T2  )/2  =  31  .3°C 
Ts  =  96. 7° C 

b.  Compute  water  velocity. 

Flow  meter  reading  =  30$  and  from  the  calibration  curve 

V  *  1.64  a/s 

W 

c.  Compute  mass  flow  rate  of  cooling  water. 

From  Reference  16,  Table  A. 6  at  304K  the  following 
properties  are  found: 
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Pr  =  5.29 

U  =  0.781  S-3  N-s/m2 
Cp  =  4.178  kJ/kg-K 

k  =  0.62  W/m-K  vw  =  1.005  S-3  m^/kg 
A  =  TtDi2/4 
A  =  1 .267  E-4  m2 
m  =  V  A/v 

W 

m  =  .2067  kg/s 

d.  Compute  heat  transferred  to  cooling  water. 

Q  =  2133  W 

e.  Compute  LMTD. 

LMTD  =  (T2-T1)/ln  [(T3-T1)/(Ts-T2)j 
LMTD  =  65.40°C 

f.  Compute  overall  heat- transfer  coefficient  based  on 

outside  area. 

A  =  -rrD  L 
o  o 

U  =  Q/A  LMTD 
o  c 

UQ  =  4086.71  W/m2-K 

g.  Compute  wall  resistance  based  on  outside  area. 

R„  ’  l-(VDi>/2  kcu 

R  =  1 .003  E-5  m2-K /W 
m 

h.  Compute  Reynolds  number  of  cooling  water. 

ae„  ■  \  “i^f 

Re  =  26535 

w 

i.  Assume  fin  efficiencies  of  =  0.20  and  n£  =  0.24 

and  assuming  (u/u,,)  ”*  ^=1  . 0 

w 
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j .  Compute  inside  heat-transfer  coefficient, 
h,  =  (k  C./D, )  Re./8  Pr*333  (u/u  J*U 


l  w  1  l  w 

=  8533.72  W/m2-K 


k.  Calculate  fin  efficiency  using  previously-calculated 
heat- transfer  coefficient  from  j  . 

=  83. 56  m  ^ 

m2  =  117.70  m"1 

n1  =  0.20 

n2  =  0.24 


1.  Calculate  cooling  water  rise. 


T  =  T,  +  AT  =  72  C 
w  b  w 


From  Reference  1 6  at  345  K,  v  =  389E-6  N-S/m‘ 

w 


(U/Uu) 


0.14 


=  1.1 


m.  Recalculate  heat  transfer  coefficient  and  iterate  if 
not  within  1  %  of  the  value  found  in  j. 

h±  =  9387  W/m2-K 

n.  Calculate  condensation  heat- transfer  coefficient. 


ho  = 


Di(L+n1L1 +n2L2^hi 


h  =  10405  W/m  -K 
o 


This  value  for  the  condensation  heat-transfer  coefficient 
agrees  within  1%  of  the  computer  generated  value  (h  =  10550) 


49 


APPENDIX  B 


SYSTEM  STARTUP  PROCEDURES 

1.  Familiarize  all  valve  numbers  by  referring  to  Figure  2.1. 

2.  Fill  boiler  with  distilled  water  to  a  level  4  inches 
below  top  of  sight  glass. 

a.  Start  vacuum  pump  or  air  ejector. 

b.  Open  suction  line  supply  valve  (3). 

c.  Allow  vacuum  pump/air  ejector  to  run  until  15  inches 
Hg  of  vacuum  is  attained. 

d.  Close  suction  line  supply  valve  (3). 

e.  Shut  down  vacuum  system. 

f.  Open  discharge  on  distilled  water  tank. 

g.  Open  distilled  water  supply  line  valve  (6). 

h.  Close  distilled  water  supply  line  valve  (6)  when 
boiler  is  filled  to  4  inches  below  the  top  of  the 
sight  glass. 

3.  For  vacuum  operations,  start  vacuum  pump/air  ejector 
then  open  suction  line  supply  valve  (3)  until  desired 
pressure  is  reached.  Suction  line  supply  valve  may  be 
throttled  to  maintain  pressure. 

4.  Energize  viewing  window  heater  and  adjust  rheostat  to 
110.  Open  air  supply  valve  (9)  about  &  turn.  The  air 
supply  may  require  further  adjustments  to  keep  viewing 
window  clear. 

5.  Open  dump  condenser  supply  valve  (2)  fully. 
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Fill  cooling  water  reservoir. 


6 . 

a.  Open  reservoir  stop  valve  (10). 

b.  Open  test  section  throttle  valve  (11). 

c.  Open  reservoir  fill  valve  and  allow  tank  to  fill  to 
about  101.6  mm  (4  in)  from  the  top. 

d.  Close  reservoir  stop  valve. 

e.  Energize  cooling  water  circulation  pump.  The  switch 
is  located  above  the  pump. 

7.  Shut  the  following  power  supply  breakers: 

a.  Panel  p-5  switch  3. 

b.  Switch  3  on  side  of  control  panel. 

c.  Switch  1  on  front  panel. 

8.  Increase  boiler  heater  voltage  to  a  maximum  of  200 
volts  during  warmup;  about  15  minutes.  When  operating  under 
vacuum  conditions  use  a  lower  voltage  to  limit  vibration 

of  glass  piping. 

9.  As  steam  is  generated,  air  is  expelled  through  a  relief 
valve.  System  pressure  should  never  be  allowed  to  exceed 
10  psig. 

10.  Once  steam  reaches  condenser  test  section  (about  30 
minutes  following  startup),  full  power  operation  may  begin. 

11.  The  apparatus  should  be  allowed  to  steam  for  30  minutes 
prior  to  taking  any  data. 
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APPENDIX  C 


SYSTEM  SHUTDOWN  PROCEDURES 

1.  Decrease  boiler  heater  voltage  to  zero. 

2.  Open  all  three  circuit  breakers. 

3.  Secure  view  window  heater  rheostat  switch;  allow  air 
to  flow  for  about  1 5  minutes  longer. 

4.  Secure  cooling  water  circulation  pump. 

5.  Close  cooling  water  throttle  valve. 

6.  Close  dump  condenser  supply  valve. 

7.  If  in  use,  secure  vacuum  suction  supply  valve  followed 
by  vacuum  pump/air  ejector. 

8.  For  cleaning  purposes  only,  the  boiler  may  be  drained 
by  unscrewing  the  Swaglok  nut  housing  the  condensate  return 
thermocouple.  The  Teflon  ferrules  require  only  hand 
tightening. 


52 


APPENDIX  D 


TEST  SECTION  INSTALLATION  PROCEDURES 

1 .  Install  outlet  nylon  holder,  mixing  chamber  and  return 
line  onto  provide  6.35  mm  (0.25  in)  diameter  bolts  and 
tighten . 

2.  Place  clean  test  tube  into  inlet  nylon  holder  and  care¬ 
fully  slide  holder  and  test  tube  into  condenser  test  section 
snugly. 

3.  Align  test  tube  to  have  one  thermocouple  channel  vertical. 

4.  Carefully  slide  thermocouples  into  channels. 

5.  Ensure  wires  are  straight  and  lined  up  with  markings  on 
inlet  nylon  holder  and  tape  into  place. 

6.  Position  nylon  thermocouple  wire  retainer  over  wires. 

7.  Place  inlet  nylon  holder  over  end  of  test  tube  and 
tighten  all  bolts. 

8.  Energize  circulation  pump  and  check  for  leakage  around 
connections  both  inside  and  outside  test  section. 


APPENDIX  E 

THERMOCOUPLE  CALIBRATION 


1 .  EQUIPMENT  USED 

a.  Thermocouple  wire 

Ccpper-constantan,  0.245  mm  (0.01  in)  Teflon  coated  wire 
was  used  for  all  thermocouples.  (Omega  Catalogue) 

b.  Calibration  bath 

A  Rosemount  Engineering  Co.  Model  9 1 3 A  calibration  bath 
shown  schematically  in  Figure  A.1  was  used. 

1)  Heating:  Electrical 

2)  Cooling:  Liquid  Nitrogen 

Note: 

Once  a  desired  temperature  is  reached,  the  temperature 
is  held  constant  by  rapid  cycling  between  heating  and  cooling. 
The  manufacturer  rates  the  temperature  bath  fluctuation 
to  be  less  than  0.01°C.  Observations  during  the  calibra¬ 
tion  procedure  showed  bath  temperature  fluctuation  to  be 
less  than  0.025°C  (I.OuV). 

c.  Thermocouple  readout 

An  HP  3054A  Automatic  Data  Acquisition/Control  System 
with  resolution  of  0.1  uV  was  used  to  obtain  data. 

d.  Reference  Point 

A  reference  junction  was  placed  in  a  0.15  m  Dia  X  0.28  m 
high  Thermos  flask  filled  with  finely-crushed  ice. 
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e. 


Temperature  measurement 


Mercury-in-glass  thermometers  with  a  resolution  of  0.1°C 
were  used.  A  total  of  four  thermometers  were  used  with  the 
following  ranges:  0-25.  25-54.  50-76  and  76-104°C. 

2.  PREPARATION 

a.  Procedure 

Two  thermocouples  were  prepared,  one  at  the  beginning 
and  the  other  at  the  end  of  a  reel  of  thermocouple  wire, 
using  the  following  procedure: 

1 )  The  Teflon  insulation  was  removed  for  a  length  of 
about  4  mm  from  one  end  of  a  2  m  long  piece  of  wire 
and  a  thermocouple  bead  was  made  using  a  Dynatech 
Corporation  thermocouple  welder. 

2)  The  other  end  of  the  thermocouple  wire  was  soft 
soldered  to  a  pair  of  heavy  duty  copper  wires. 

3)  Heat  shrink  insulation  was  placed  on  the  soldered 
end. 

b.  Reference  point 

The  thermocouple  reference  junctions  were  placed  in  a 
5  cm  test  tube  and  then  into  a  thermos  flask  containing 
finely  crushed  ice  as  shown  in  Figure  A.1. 

c.  Thermometers 

Thermometers  were  completely  immersed  in  the  calibration 
bath  as  recommended  by  the  manufacturer. 

d.  Analysis 

A  short  computer  program  (TC_CAL)  was  written  to  accept 
the  thermocouple  readings  through  the  HP  data  acquisition 
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system,  and  the  bath  temperature  through  the  keyboard.  This 
program  prints  all  data  as  well  as  stores  the  data  on  a 
computer  disk.  A  listing  of  the  computer  program  is  provided 
in  Appendix  G. 

3.  CALIBRATION  PROCEDURE 

a.  The  bath  temperature  was  set  at  about  15°C. 

b.  A  constant  temperature  is  reached  when  rapid  cycling 
of  heating  and  cooling  occurs,  but  wait  three  minutes  to 
enter  data. 

c.  Enter  the  bath  temperatures,  as  measured  by  the 
mercury-in-glass  thermometer,  into  the  keyboard.  Immediately 
following  the  data  entry,  the  computer  will  read  the 
thermocouple  EMF  automatically  and  print  the  bath  temperature 
and  thermocouple  EMF  (yV). 

d.  When  the  computer  prompts  for  more  data,  increase 
bath  temperature  by  about  10°C  and  repeat  the  above  procedure 
up  to  about  95°C. 

e.  Repeat  the  above  procedure  for  decreasing  temperatures. 

4.  CALIBRATION  CURVE 

a.  A  least-squares  technique  was  used  to  generate  a 
third  order  polynomial  in  the  form  of: 

T  =  a^ E+a2S^+a^S^ 

T  is  the  temperature  in  degrees  Celsius  and  E  is  the 
thermocouple  reading  in  microvolts. 

b.  Values  generated  for  the  coefficients  are: 
a1  =  0.02635 
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a.j  =  6.  577  E-11 

c.  Figure  A. 2  is  a  sample  printout  of  calibration 
data. 

d.  Figure  A. 3  is  a  graphical  representation  of  calibra¬ 
tion  data. 


Figure  E.1  Schematic  of  Calibration 


Month,  date  and  t  ime  :  03  : 1  6 : 20  : 20  : 30 


Sat'n  Temp 

Enf (0> 

En+  <  1  > 

(Deg  C> 

(Micro  V) 

(Micro  V> 

13.06 

503.0 

503.2 

24.33 

955.0 

956.0 

35 . 33 

1409.5 

1409.0 

46.72 

1889.8 

1891 . 1 

59.50 

2442.9 

2442.4 

72.06 

2996.4 

2995.7 

34.78 

3575.8 

3575.  ! 

94.00 

3993.8 

3998 . 1 

93.22 

3963.9 

3963.0 

82.89 

3483.8 

3488.2 

66.89 

2767.8 

2766.9 

54.33 

2216.9 

2215.8 

42.06 

1691 .7 

1690.9 

30.67 

1213.4 

1212.4 

20.72 

806.5 

806.0 

13.06 

502.5 

502.0 

16  runs  were  stored  in  file  CAL_DflTA 


Figure  E.2  Calibration  Printout. 
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ITS 


Thermocouple  Calibration  Curve. 


APPENDIX  F 


TEST  TUBE  CLEANING  PROCEDURES 
- Caution - 

The  following  mixture  is  highly  toxic  and  very  irritant 
when  contacted  by  the  skin.  Be  very  careful  not  to  splash 
the  solution  in  any  way.  Wear  protective  gloves  and  eye 
protection  when  the  actual  cleaning  is  being  conducted. 

1 .  Ensure  the  surface  of  the  tube  is  smooth  by  placing 
in  lathe  and  polishing  with  copper  cleaner  (Brasso)  and  a 
clean  rag. 

2.  While  still  in  the  lathe,  wipe  completely  with 
acetone  followed  by  ethyl  alcohol. 

3.  Prepare  a  solution  of  equal  parts  of  5056  sodium 
hydroxide  solution  (commercially  prepared,  available  through 
MCB  Manufacturing  Chemists  Inc.,  Cincinnati,  OH)  and 

ethyl  alcohol.  About  4  ounces  each  mixed  in  a  shallow 
stainless-steel  pan  is  sufficient.  The  sodium  hydroxide 
will  precipitate  out  of  solution  and  form  a  white  paste. 

4.  Heat  the  mixture  to  about  80°C,  but  try  to  avoid 
spattering.  The  white  pasty  texture  os  the  solution  should 
now  dissolve  into  a  clear  liquid. 

5.  Completely  immerse  the  tube  into  the  mixture  and 
scrub  the  entire  surface  with  a  bristle  brush  (old  tooth¬ 
brushes  work  fine). 

6.  Rinse  with  tap  water  followed  by  distilled  water. 
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7.  Allow  the  tube  to  remain  in  distilled  water  until 
ready  for  installation. 

8.  The  nylon  holders  should  also  be  cleaned  with  the 
mixture . 

9.  The  test  tube  should  be  installed  as  soon  as  possible 
after  cleaning. 


APPENDIX  G 


COMPUTER  PROGRAM  LISTINGS 


100!  FILE  NAME:  TC_CAL 
110!  REVISED:  June  1 .  1983 

120! 

130  DIM  En.f <  1  > 

140  PRINTER  IS  701 
ISO  BEEP 

ISO  INPUT  '‘ENTER  MONTH.  DATE  AND  TIME  ( MM : DD : HH : MM : SS  > " . DateS 
170  OUTPUT  709 : "TD“ :DateS 
ISO  OUTPUT  709 : "TD" 

190  ENTER  709: DateS 

200  PRINT  USING  "12X .""Month.  date  and  t i me 1 4A“ : DateS 
210  BEEP 

220  INPUT  "GIVE  A  NAME  FOR  DATA  FILE" ,D_f 1 leS 

230  CREATE  BDAT  D_fileS.5 

240  ASSIGN  »File  TO  D_fileS 

250  J-0 

2S0  PRINT  "  " 

270  PRINT  USING  "l2X.""Bath  Temp  Ernf(O)  Emfd) . 

280  PRINT  USING  "12X.“"  (Deg  C)  (Micro  V)  (Micro  V) . 

230  Repeat:  ! 

300  BEEP 

310  INPUT  "ENTER  BATH  TEMPERATURE  (DEG  F)",T_bath 
320  T_bath-(T_bath-32. )/1 .8 
330  OUTPUT  709 ; "AR  AF20  AL21" 

340  OUTPUT  722: "FI  R1  T1  Z1  FL1 " 

350  FOR  1-0  TO  1 
3SC  OUTPUT  709: "AS  SA" 

370  ENTER  722:Emf(I> 

380  Emf (I)-ABS(Emf (I))»10“6 
390  NEXT  I 

400  PRINT  USING  " 1 3X . 3D . DD . 7X , 4D . D . 5X . 4D . D" : T_bath .Emf (*) 

410  OUTPUT  “File;  T_bath ,Emt‘ <«) 

420  3EEP 

430  INPUT  "ARE  YOU  TAKING  MORE  DATA  ( 1  - YES . O-NO) ?" . Go_on 
440  J-Jvl 

450  IF  Go_on* 1  THEN  Repeat 

4bU  BttP 

470  PRINT  "  " 

480  PRINT  USING  "I2X.DD.""  runs  were  stored  In  file  "”  .  I OA"  :  J ,D__f  i  leS 
490  END 


’000!  FILE  NAME:  NSFDRP 
1010  ! REVISED:  23  MAY  1983 
’020  DIM  Emf(IO) 

’030  Di* . 01 27 
'040  Do*. 01905 
1050  Dr*. 015875 
1060  Dssp- . 1 524 

1070  Ax-Pl*Dssp '2/4-PI»Do»L  !  TO  EE  MODIFIED 

1080  L-. 13335 

’090  LI -.060325 

1100  L2*. 034925 

1110  L 1 -L 1 +<  Do-Di ) /4 

1120  L2*L2+(Dr-Di)/4 

1130  Ci*. 029 

1140  Kcu-385 . 

1150  Rm»Do*LOG< Do/D i )/<2»Kcu) 

1160  PRINTER  IS  701 
1170  CLEAR  709 
1180  BEEP 

M90  INPUT  "ENTER  MONTH.  DATE  AND  TIME  <MM:DD:HH:MM:SS>" .DateS 
1200  OUTPUT  709:"TD" ;Date$ 

1210  OUTPUT  709: "TD" 

1220  ENTER  709;Date$ 

1230  PRINT  "  Month,  date  and  tine  :":DateS 

1240  BEEP 

1250  INPUT  "ENTER  INPUT  MODE  < 1 -3054A .2-FILE)" . Ira 
1260  IF  Im-1  THEN 
1270  BEEP 

1280  INPUT  "GIVE  A  NAME  FOR  THE  RAW  DATA  FILE"  ,D_f i ieS 
1290  CREATE  BDAT  D_fileS.10 
1300  ELSE 
1310  BEEP 

1320  INPUT  "GIVE  THE  NAME  OF  THE  EXISTING  DATA  FILE" . D_f i IeS 
1330  PRINT  "  " 

1340  PRINT  USING  "12X.""The  following  analysis  was  performed  for  data  in  file 
". 10A":D_f ileS 
’350  BEEP 

1360  INPUT  "ENTER  THE  NUMBER  OF  RUNS  STORED" .Nrun 
1370  END  IF 
1380  BEEP 

1390  INPUT  "GIVE  A  NAME  FOR  PLOT  DATA  FILE" .P_f i IeS 
1400  CREATE  BDAT  P_fileS.5 
1410  ASSIGN  *File  TO  D_fileS 
1420  ASSIGN  SFilep  TO  P_file$ 

1430  J-0 

1440  Repeat:! 

1450  J-J+1 

1460  IF  Ira-1  THEN 

1470  BEEP 

1480  INPUT  "ENTER  FLOWMETER  READING". Fra 
1490  OUTPUT  709:"AR  AF20  AL30" 

1500  OUTPUT  722:"F1  R1  T1  Z1  FU" 

1510  OUTPUT  709:"AS  SA" 

1520  ENTER  722:Bpv 
1530  FOR  1-0  TO  10 
1540  OUTPUT  709: "AS  SA" 

1550  IF  I >4  THEN 

1560  Se-0 

1570  FOR  K-0  TO  10 
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1580  ENTER  722 ;E 
1590  Se-Se+E 
>600  NEXT  K 

1610  Emf  <I)*A8S(Se/10)»1 .£+6 

1620  ELSE 

1 630  ENTER  722  ;E 

1640  Emf(I)-A8S(E>»1 .£+6 

1650  END  IF 

1660  NEXT  I 

1670  ENTER  722:Emf(I> 

1680  OUTPUT  71 3:''T1R2E" 

1690  WAIT  2 

1700  ENTER  713;T1 

1710  OUTPUT  7 13:"T2R2£" 

1720  WAIT  2 
1730  ENTER  713; T2 
1740  IF  J-1  THEN 
1750  BEEP 

1760  INPUT  "ENTER  NANOMETER  READING  (mm  Hg)",Phg 
1770  BEEP 

1730  INPUT  "ENTER  HEIGHT  OF  MANOMETER  WATER  COLUMN" .Pwater 
1790  END  IF 
1800  ELSE 

1310  ENTER  9FiLe:Emf<*> .Fm.TI .T2.Phg. Pwater. Bpv 
1820  END  IF 
1830  Tsteam-0 
1840  FOR  1-0  TO  1 

1 850  Tsteam-Tsteam+ . 5*FNTvsv( Emf ( I >  > 

1860  NEXT  I 
18 70  Twm-O. 

1830  FOR  1-0  TO  5 

1890  Tut  I ) -FNTvsv(Emf  < 1+5) > 

1900  Turn* Twm+Tw( I > 

1910  NEXT  I 

1920  Twm-Twm/6 

1930  Tcon-FHTvsv<Emf < 10) ) 

1940  Psat-FNPvst ( Tsteam) 

1950  P  tes  t - < Phg+Puate r / 1 3 . 6  >  *  1 33 . 322 

1960  IF  J-1  THEN 

1970  Vst-FNVvst< Tsteam) 

1980  Ppng*(Ptest-Psat)/Ptest 
1990  Ppst* 1 -Ppng 

2000  Mt  ng-Ppng/(287»<  Tsteam’-273 .  15) / (Vst  "P  test  >  +Ppng) 

2010  Of ng-Mf ng/( 1 . 608- . 608»Mf  ng) 

2020  Mfng-Mfng*100 
2030  Of  ng-Vf  ng*1 00 
2040  BEEP 
2050  PRINT  "  " 

2060  PRINT  USING  "  1 2X , ’"'Measured  pressure  * 

E.""  <Pa) . :Ptest 

2070  PRINT  USING  *’  1 2 X , ""Pressure  corresponding  to  measured  steam  temp  * 

E.  ""  (Pa) . :Psat 

2080  °RINT  USING  "12X, ""Noncondensable  gas  concentration  -"".D.l 

Vo  I  '/.) . :  Vf  ng 

2090  PRINT  USING  "12X. ""Noncondensable  gas  concentration  •"" .DD 

Mass  X) . :Mfng 

21 00  IF  Mf ng> . 5  THEN 
2110  BEEP 
2120  PRINT  "  " 

2130  IF  In-1  THEN 


.D.3D 
,D.3D 
.""  ( 
D .""( 
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2140  BEEP 
2150  PRINT  "  " 

2160  PRINT  USING  ”1  OX , ’"'Energi le  the  vacuum  system  . 

2170  BEEP 

2130  INPUT  "OK  to  accept  the  present  run  < t -YES ,0*N0)^" . Ok 
2190  IF  Ok-O  THEN 
2200  8EEP 

2210  DISP  "NOTE:  THE  PRESENT  DATA  RUN  WILL  BE  DISCARDED!!  " 

2220  WAIT  5 

2230  GOTO  1460 

2240  END  IF 

2250  END  IF  ' 

2260  END  IF 
2270  END  IF 

2280  IF  Im-t  THEN  OUTPUT  3Fi IeiEraf (»> ,Fn. T ! . T2.Phq,Pmater ,Bos 
2290!  ANALYSIS  BEGINS 
2300  Tavg*  <  T 1 +T2) *  .5 
2310  Cpu*FNCpu< Tavg) 

2320  Rho<J*FNRhou<  Tavg) 

2330  Mf-Fm*. 02451*. 3043*3/100 

2340  Md*f1f*Rhou 

2350  Vu-Mf/<PI»Di"2/4) 

2360  Q-Md»Cpui»<  T2-T1  > 

2370  Op-O/ <  PI»Do»L) 

2380  Two-Tum+Qp-Rm-  .5 
2390  Dtf-Tsteam-Tuo 
2400  Ku-FNKu< Tavg) 

2410  Muu*FNMuui<  Tavg) 

2420  Reu-Rhou-Vu-Di/Muw 
2430  Pru-FNPru< Tavg) 

2440  Fe 1 *0 . 

2450  Fe2-0. 

2460  Cf-1. 

2470  Hi *Kij»Ci /D  i -Reu  ' . 3-Pru  ' .  3333 -Cf 
2480  Dt-Q/<PI-Di-<L+L1-Fe1+L2-Fe2>-Hi> 

2490  Cf  c* <Muu/FNMuu< Tavg+D t  >>".14 

2500  IF  ABS<  <Cfc-Cf)/Cfc>>.01  THEN 

2510  Cf-(Cf+Cfc)-.5 

2520  GOTO  2470 

2530  END  IF 

2540  PI -PI-(Dl+Do) 

2550  A!-<Do-Di>-PI-<Di+Do>».5 

2560  Ml  * < H i -PI / < Kcu-A1 ) ) A .5 
2570  P2-PI-<Di+Dr> 

2580  A2-<Dr-Di)»PI»<Di+Dr)".5 

2590  M2”<Hi-P2/<Kcu-A2> ) " .5 

2600  Lmtd*<  T2-T1 ) /L0G< < Tsteam-TI > /< Tsteam-T2> > 

2610  Uo-G/<Lmtd»PI«Do»L) 

2620  Fel *FNTanh<M1«L1 >/<M1-L1 ) 

2630  Fe2-FNTanh<M2-L2)/<H2-L2> 

2640  Ho-1  /<  1  /Uo-Do»t/<Dl-<L+Li  -Fel  +12-Fe2)-Hi>-Rm> 

2650  Dtc-0/<PI-Di-(L+L1-Fe1+L2-Fe2>-Hi> 

2660  IF  ABS<<Dtc-Dt>/Dtc)>.01  THEN  2470 
2670  T room*FNTvsv< (Emf  <3) +Emf  <  4) ) -.5) 

2680  PRINT  "U0-" :Uo 

2690  Hos*Q/<°I-Po-L-< Tsteam-Two) ) 

2700  PRINT  "HO  STAR-";Hos 
2710  Hf  g*FMHf  g<  Ts team) 

2720  Tf t lm*< Tsteam+Tuo)-.5 
2730  Kf'FNKu<Tf ilm) 
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2740 

2750 

2760 

2770 

2780 

2780 

2800 

2810 

2820 

2830 

2840 

2850 


Rhof  -FNRhoui T  film) 
Muf  -FNMuw<  Tf  i  lm> 


Hnu-.725*<Kf  '3-Rhof '2-9.81»Hfg/<Huf-Do*Dtf )>' .25 
Hpq-  .  65 1  »Kf  «<  Rhof  *2»9 . 81  «Hf  g/(Muf  *Do*0p>  > '  .3333 
PRINT  M  " 

PRINT  USING  "12X. ""Water  velocity  -  "" 

PRINT  USING  "12X. ""Reynolds  number  -  "" 

PRINT  USING  " 12X Inside  heat-tran  coef 
PRINT  USING  "1 2X .""Heat  flux  -  "" 

PRINT  USING  "12X. ""Outside  heat-tran  coef  *  "" 

PRINT  USING  "l2X.""Nusselt  coefficient  -  "" 

PRINT  USING  "12X. 12X,""Nus  coef  (from  heat  flux) 


Z .  DD . ""  (m/S> . ;Vu 

5D.D":Reu 

5D.D.""  (N/n“2-K> . ;Hi 

Z.5DE."”  <W/n'2> . ;0p 

5D . D . ""  <W/m  2-K> . :Ho 

5D . D . ""  <W/n'2-K> . ;Hnu 

*  ""  ,5D>D ,  ""<W/r»  "2-K> . 


Hqp 

2860  OUTPUT  »Fi lep :0p . Hos .Hnu 
2870  IF  Im-1  THEN 
2880  BEEP 

2890  INPUT  "ARE  YOU  TAKING  MORE  DATA  ( 1 -YES.O-NO) ?" ,Go_on 
2900  IF  Go_on-1  THEN  Repeat 
2910  ELSE 

2920  IF  J<Mrun  THEN  Repeat 
2930  END  IF 
2940  IF  Im-1  THEM 
2950  BEEP 

2S60  PRINT  USING  "J2X,""N0TE:  "" ,DD , ""  data  runs  were  stored  in  file  "".10A":J. 

D_f ileS 

2970  END  IF 

2980  BEEP 

2390  PRINT  "  " 

3000  PRINT  USING  "12X.""N0TE:  "" .DD .""X-Y  pairs  were  stored  in  plot  data  file  " 

". 10A";J.P  files 

3010  ASSIGN  9F  i  1 e  TO  • 

3020  END 

3030  DEF  FNPvst < Tsteam) 


3040  DIM  K < 8 > 

3050-  DATA  -7.691 234564 . -26 . 08023696 .-158.1 706546 . 64 . 23285504 .-113. 9646225 
3060  DATA  4 . 1671 1 732 ,20 .9750676 , 1 .E9 ,6 
3070  READ  K<»> 

3080  T-<  Ts team +273. 15>/647 .3 

3090  Sum-0 

3100  FOR  N-0  TO  4 

3110  Su/n-Sum+K(N>*(  I  -  D  “(N+1 ) 

3120  NEXT  N 


3130 


3r-Sum/<  T-( 1 +K(5>-<  1  -T>  i 


K(G>*( 1-T>  2))-<!-T)/<K(7>*('-T) '2+K<8)> 


3140  Pr-EXP(Br) 

3150  P-22l20000»Pr 

3160  RETURN  P 
3170  FNEND 
3180  DEF  FNHfg<  T> 

3190  Hfq-2497 . 7389-T»<2. 2074+T»< 1 . 7079E-3-2 . 8593E-6) > 

3200  RETURN  Hfg-1000 

3210  FNEND 

3220  DEF  FNMuu(T) 

3230  Mu-1 . 5760947 3E-3-TM3. 51 198576E-5-T-<3. 58358 1 6E- 7- 1 . 365586 1 1 5E-9»T  > ) 

3240  RETURN  Mu 

3250  FNEND 

3260  DEF  FNVvst(T) 

3270  V-58 . 4525538-T«(  1 .5 1 508776-W  .  01 37274658S-T-4 .253667 1 1E-5> ) 

3280  RETURN  V 
3290  FNEND 
3300  DEF  FNCpufT) 
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3310  Cpw-4.21 120858-T*<2.26826E-3-T*<4.42361E-5+2.71428E-7«T> ) 

3320  RETURN  CpuMOOO 

3330  FNEND 

3340  DEF  FNRhou< T) 

3350  Ro*999 . 52946 +T«< .01269-T*<5.482513E-3-T*1 .234147E-5)) 

3360  RETURN  Ro 
3370  FNEND 

3380  DEF  FNPrw< T )  „ 

3390  Prw-IO'O . 09976605-T«< 1 .3749326E-2-T*<3.9S8875E-5-3.45026E-7»T>  >> 

3400  RETURN  Prw 

3410  FNEND 

3420  DEF  FNKu<T) 

3430  Ku-.5625894+T«(2.2964546E-3-T«< 1 . 509766E-5-4 . 0581652E-8-T) > 

3440  RETURN  Kw 

3450  FNEND 

3460  DEF  FNTanh(X) 

3470  P-EXP(X) 

3480  Q-EXP(-X> 

3490  Tanh-(P+Q)/<P-Q) 

3500  RETURN  Tanh 
35  I  0  FNEND 

3520  DEF  FNTvsv< V >  ,  .  „ 

3530  T-VM . 0261733441 6- V»<9 .244 7859E- 7 -U»6. 0746b42E~ 1 1 ) > 

3540  RETURN  T 
3550  FNEND 
3560  DEF  FNHf <T) 

3570  Hf-T*(4.203849-T«<5. 881 32E-4-T«4 . 5516031 7 E-6> ) 

3580  RETURN  Hf«1000 
3590  FNEND 


1000!  FILE  NAME :  SIEDER 
1010!  REMISED:  May  23.  1983 
1020  DIM  Entf(IO)  ,Tu)<5) 

1030  Kcu-385 

1040!  Assign  Geometric  Variables 
1050! 

1060  Di-.0127  !Inside  Diameter 

1070  Do”. 01905  fOutside  Diameter 
1080  Dr”. 015875  !0u.tlet  End  Diameter 
1090  L”. 13335  !Condensing  Length 

1100  LI”. 060325  !Inlet  End  "fin”  Length 
1110  L2-. 034925  JOutlet  End  "fin”  Length 

1120  PRINTER  IS  701 
1130  BEEP 
1140  CLEAR  709 

1150  INPUT  "ENTER  MONTH,  DATE  AND  TIME  <MM:DD:HH:MM:SS>",B$ 

1160  OUTPUT  709 ;"TD” ;BS 
1170  Series'! 

1180  OUTPUT  709;"TD" 

1190  ENTER  709: AS 

1200  PRINT  USING  "  1 2X .  "".Month .  date  and  time:  "".14A":AS 
1210  BEEP 

1220  INPUT  "ENTER  INPUT  MODE  < 1 -3054A , 2-FILE >‘\Im 
1230  IF  In-1  THEN 
1240  BEEP 

1250  INPUT  "PROVIDE  NAME  FOR  THE  DATA  F ILE" . D_f i leS 
1260  CREATE  BOAT  D_fileS,10 
1270  ELSE 
1280  BEEP 

1290  INPUT  "PROVIDE  NAME  OF  THE  DATA  FILE" ,D_f i le$ 

1300  BEEP 

1310  INPUT  "ENTER  THE  NUMBER  OF  RUMS  STORED". Nrun  r 

1320  PRINT  USING  "12X.""The  following  analysis  was  performed  for  data  in  file 
". 10A":D  files 
1330  END  IF 
1340  BEEP 

1350  INPUT  "PROVIDE  NAME  FOR  PLOT  DATA  FILE". Plots 

1360  CREATE  BDAT  Plots, 5 

1370  ASSIGN  SFile  TO  D_fileS 

1380  ASSIGN  sFilep  TO  Plots 

1390  J-0 

i/ifin  ^vs»o 

1410  Sxy-0 

1420  IF  In-1  THEN 

1430!  READ  DATA  THROUGH  DATA  ACQUISITION  SYSTEM 
1440!  IF  THE  INPUT  MODE  (In)  ”  1 
1450  BEEP 

1460  INPUT  "ENTER  FLOWMETER  READING" .Fm 
1470  OUTPUT  709: "AR  AF20  AL30" 

1430  OUTPUT  722 ; "F 1  R1  T1  Z1  FL1" 

1490  FOR  1-0  TO  10 
1500  OUTPUT  709; "AS  SA" 

1510  IF  I>4  THEN 
1520  Se-0 
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1 


1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1630 

1690 

1700! 

1710! 

1720 

1730 

1740! 

1750 

1760! 

1770 

1780 

1790 

1800 

1810 

1820 

1330! 

1840 

1850 

1860 

1870 

1880 

'890 

1900 

1910 

1920 

1930 


FOR  K-0  TO  9 
ENTER  722 :E 
Se*Se+E 
NEXT  K 

Emf  < I)*ABS(Se/1 0) 

ELSE 

ENTER  722  :E 
Emf  < I ) -ABS(E) 

END  IF 
NEXT  I 

OUTPUT  713;"T1R2E" 

WAIT  2 
ENTER  7 1 3 ;  T 1 
OUTPUT  7  T 3 ; ** T2R2E** 

WAIT  2 

ENTER  713 :T2 
ELSE 

READ  DATA  FROM  USER-SPECIFIED  FILE  IF 
INPUT  MODE  dm)  *  2 
ENTER  SFileiEmf(-) ,T1 ,T2.Fro 
END  IF 

Compute  Average  Bulk  Water  Temperature 
Tavg*  <  T1 +T2 ) *  .  5 

Compute  Average  Wall  Temperature 

Tual 1*0 

FOR  1-5  TO  10 

Tu<  1-5 ) *FNTvsv<  Emf < I) ) 

Tual  l-Tuall  +  Tud-5) 

NEXT  I 

Tual l*Twal 1/6 

Compute  The rmophys leal  Properties 
Cpu*FNCpu<  Tavg) 

Rhou-FNRhou( Tavg) 

Ku-FNKu< Tavg)  v 

Hf=.Fm-.02451».3048"3/100  !Volume  t  lourate  of  Water  <m  i/s) 

Md*Mf »Rhou  ! Mass  Flourate  of  Water  <kg/s> 

Vu»Mf/(PI*Di ’2/4)  ! Water  Velocity  (m/s) 

0-Md*Cou»( T2-T1 )  !Heat  Transferred  to  Cooling  Water  (W> 
Dtu-0*L0G(Do/Di )/<2*PI*Kcu»L)  !Temp  Drop  Across  Tube  Wall  (Deg  u) 
Tual l* Tual 1 -Dtu» .5 

Lmrd-(T2-T1 >/L0G((Tuall-T! )/<Tuall-T2>> 

P1*PI*2»Do  !Heat-Transfer  Perimeter  At  Inlet 
P2-PI»2*Dr  !Heat-Transfer  Perimeter  At  Exit 


•  ( Do-D i ) *P I*Do  ! Heat 


Area  At  Inlet 


A2*<Dr-Di )*PI»Dr  !Heat-Transfer  Area  At  Exit 
Hi-0/(PI»Di<L»Lmtd) 

Ml  * (H i»P1 / (Kcu»A 1 ) ) ' .5 

M2-(Hi-P2/<Kcu*A2)>‘.5  npr  _  T  .  . 

Fe1-FNTanh(M1*L1 )/<M1*L1 >  !Fin  Efficiency  At  Inlet 

Fe2*FNTanh(M2»L2)/<M2*L2>  !Fin  Efficiency  At  Exit 

Hic-Q/(PI*Di»(L+L1 *Fe1 +L2*Fe2)»Lmtd> 

IF  ABS<(Hi-Htc)/Hic)>.01  THEN 

Hi  - <Hlc+Hi )» .5 

GOTO  1990 

END  IF 

PRINT  "  " 

PRINT  USING  ” l2X.""Posi tlon  number  :  T 


PRINT  USING 
PRINT  USING 


•  l2X.""Posi tlon  number 

'12X.""Wall  temperature  (Deg  C) 
'12Xi’"’Average  wail  temperature 


M",6(DD.DD. IX)” :T«(*> 
””  , DD . DD (Deg  C)"’’ 


I 


2120  PRINT  USING  "12X . '"'Inlet  water  temperature  -  "" .DO . DD . 

2130  PRINT  USING  "1 2X , ""Out let  water  temperature  «'  "" . DD . DD . "" 

2140  PRINT  USING  " 1 2X . ""Log-mean- temp  difference  -  "" ,DD.DD. "" 

2150!  CALCULATE  THE  NUSSELT  NUMBER 

2160  Nu-H ic»Di /Kw 

2170  Re-Rhow-Uw-Di/FNMufTavg) 

2180  Cf - ( FNMu<  Tavg ) /FNMu< Twal 1 ) > - . 1 4 
2190  Prw-FNPrw< Tavg) 

2200  X-Re ' . 8»Prw " . 3333»Cf  ICompute  Sieder-Tate  Parameter 
2210!  COMPUTE  COEFFICIENTS  FOR  THE  LEAST-SQUARES-FIT 
2220!  STRAIGHT  LINE 

2230  PRINT  USING  "1 2X . '"’Water  velocity  -  "".Z.DD."" 

2240  PRINT  USING  "12X. '"’Sieder-Tate  parameter  -  "",5D.D":X 

2250  PRINT  USING  "12X.""Nusselt  number  -  "".4D.DD";Mu 

2260  OUTPUT  <*Filep:X,Nu 

2270  Sx»Sx+X 

2280  Sy-$y *Nu 

2290  Sxs-Sxs+X<*X 

2300  Sxy-Sxy+X*Nu 

2310!  STORE  RAW  DATA  IN  A  USER-SPECIFIED  FILE  IF 
2320!  INPUT  MODE  (In)  »  1 

2330  IF  In-1  THEN  OUTPUT  9File;Emf<«> ,T1 ,T2,Fm 

2340  BEEP 

2350  J-J+1 

2360  IF  In-1  THEN 

2370  INPUT  "ARE  YOU  TAKING  MORE  DATA  < 1 -YES , 0-NG>?" ,Go_on 
2380  Nrun-J 


(Deg  C) . : Tt 

(Deg  C>"""i T2 
(Deg  C) . ;Lmtd 


(m/S> . ;Vw 


2390  IF  Go_on- 1  THEN  1420 
2400  ELSE 

2410  IF  J<Nrun  THEN  1420 
2420  END  IF 
2430  Cx-Sxy/Sxs 
2440  PRINT  "  " 

2450  PRINT  USING  "12X. ""Sieder-Tate  Coefficient  -  "",D.4D":Ci 
2460  Ac-0. 

2470  PRINT  "  " 


:4.90  PRINT  USING  "12X. ""Least-Squares  Line: . 

""Slope  •  ”".MD.SDE.":Cl 


2490 

PRINT  USING  "  1  4X , 

2500 

PRINT  USING  " 1 4X , ' 

2510 

IF  Im-1  THEN 

2520 

BEEP 

2530 

PRINT  USING  "12X2 

n.D  f 

ileS 

2540 

run  rr 

l_l  It/  X  I 

2550 

ASSIGN  @File  TO  * 

2560 

ASSIGN  @F i lep  TO 

2570 

END 

2530 

DEF  FNRhow(T) 

2590 

Ro- 1 006 .35724- T»< 

2600 

RETURN  Ro 

2610 

FHEND 

2620 

DEF  FNPrw(T) 

2630 

Pr-10  *  < 1 .09976605 

2640 

RETURN  Pr 

2650 

FNEND 

2660 

DEF  FNMu<  T) 

2670 

Mu-1 .57609473E-3- 

2630 

RETURN  Mu 

2690 

FNEND 

2700 

DEF  FNKw< T) 

1  .DD. 


"" .MD.5DE ! Ac 

data  runs  were  stored  in  file  "”.3A":Nru 
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2710  Kw* .572133504477+1 . 52770121209E-3*T 

2720  RETURN  Kw 

2730  FNEND 

2740  DEF  FNTvsv(Enf) 

2750  V-Emf*10'6 

2760  T-V-< .026 173344 16-V**9. 2447859E-7-V*6. 0746642E-1 1>> 

2770  RETURN  T 
2780  FNEND 

2800  Cpw-<4?2U20858-T»<2 . 26326E-3-T*<4 . 4236  IE-5 +  T»2 . 71 428E-7) ) >*1000 

2810  RETURN  Cpu 

2820  FNEND 

2830  DEF  FNTanMX) 

2840  P-EXP<X> 

2850  O-EXP(-X) 

2860  Tanh-(P+Q)/<P-Q> 

2870  RETURN  Tanh 
2880  FNEND 
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(000!  FILE  NAME:  PLOT 
1010  PRINTER  IS  705 
1020  BEEP 

1030  INPUT  "ENTER  MINIMUM  AND  MAXIMUM  X-VALUES" .Xmln .Xmax 
1040  BEEP 

1050  INPUT  "ENTER  MINIMUM  AND  MAXIMUM  Y-VALUES" . Ymin , Ymax 
1060  BEEP 

(070  INPUT  "ENTER  STEP  SIZE  FOR  X-AXIS" .Xstep 
1080  BEEP 

1080  INPUT  "ENTER  STEP  SIZE  FOR  Y-AXIS" .Ystep 
1100  BEEP 

TITO  PRINT  "INiSPI : IP  2000.2000.8000.7000;" 

1120  PRINT  "SC  0. 100.0. 100;TL  2.0;" 

(130  Sf X- ( 00/ ( Xnax-Xm in > 

1140  Sfy-100/( Ymax-Ym in) 

1150  PRINT  "P'J  0.0  PD" 

1160  FOR  Xa-Xmin  TO  Xmax  STEP  Xstep 

(170  X-<Xa-Xmin)*$fx 

1180  PRINT  "PA" ; X . " . 0 ;  XT;" 

(190  NEXT  Xa 

1200  PRINT  "PA  100.0:PU:" 

(210  PRINT  "PU  PA  0.0  PD" 

1220  FOR  Ya-Ymin  TO  Ymax  STEP  Ystep 
(230  Y-<Ya-Ymin)*Sfy 
1240  PRINT  "PA  0,";Y,"YT" 

1250  NEXT  Ya 

1260  PRINT  "PA  0.(00  TL  0  2" 

1270  FOR  Xa-Xmin  TO  Xmax  STEP  Xstep 

(280  X-<Xa-Xmin)*Sfx 

1290  PRINT  "PA":X,",100;  XT" 

1300  NEXT  Xa 

1310  PRINT  "PA  100.100  P'J  PA  100.0  PD" 

1320  FOR  Ya-Ymin  TO  Ymax  STEP  Ystep 

1330  Y-<Ya-Ymin)-Sfy 

1340  PRINT  "PD  PA  100,",Y,"YT" 

1350  NEXT  Ya 

1360  PRINT  "PA  100.(00  PU" 

1370  PRINT  "PA  0,-2  SR  1 .5.2" 

1380  FOR  Xa-Xmin  TO  Xmax  STEP  Xstep 
1390  X-(Xa-Xnin)*Sfx 
! 400  PRINT  "PA":X,".0:" 

1410  PRINT  "CP  -2,-1 ;LB";Xa;"" 

1420  NEXT  Xa 

1430  PRINT  "PU  PA  0.0" 

1440  FOR  Ya-Ymin  TO  Ymax  STEP  Ystep 
1450  Y-< Ya-Ymln)»Sfy 

(460  PRINT  "PA  0.":r."" 

1470  PRINT  "CP  -4.-.25;LB";Ya;"" 

1480  NEXT  Ya 
1 490  BEEP 

1500  INPUT  "ENTER  X-LABEL" ,X labels 
1510  BEEP 

1520  INPUT  "ENTER  Y-LABEL" .Y labels 

1530  PRINT  "SR  1.5.2;PU  PA  50.- 10  CP";-LEN<XlabelS)/2;"0:LB":XIabelS:"" 
1540  PRINT  "PA  -11.50  CP  0 , " : -LEN< Y labe IS > /2*5/6 ; "DI  0 . 1 :LB" ; YlabelS;"" 
1550  PRINT  "CP  0.0" 

1560  Repeat:! 

1570  BEEP 
V 
y 


1580 
1590 
1S00 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
21  10 
2120 
2130 


INPUT  "ENTER  THE  NAME  OF  THE  DATA  FILE".D  fileS 
ASSIGN  $File  TO  D_file$ 

BEEP 

INPUT  "ENTER  THE  BEGINNING  RUN  NUMBER". Md 
BEEP 

INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED" .Npairs 
BEEP 

INPUT  "SELECT  A  SYMBOL  FOR  THE  PLOTTER  <1-*,2-+.3-c,4-o.5-')",Sy 
PRINT  "PU  DI" 

IF  Sy-1  THEM  PRINT  "SM*" 

IF  Sy-2  THEN  PRINT  "SM+" 

IF  Sy-3  THEN  PRINT  "SMc" 

IF  Sy-4  THEN  PRINT  "SMo" 

IF  Sy-5  THEN  PRINT  "SM“" 

IF  Md>Npairs  THEN 
FOR  1-1  TO  (Md- 1 ) 

ENTER  ftFile;Xa.Ya 
NEXT  I 
END  IF 

FOR  1-1  TO  Npairs 
ENTER  <3Fiie;Xa.Ya 
X-(Xa-Xmin)«Sfx 
Y-<Ya-Ymin)»Sfy 
PRINT  "PA"  ,X,  Y,  '’* 

NEXT  I 
BEEP 


ASSIGN  3F,le  TO  - 

INPUT  "DO  YOU  HAVE  MORE  DATA  TO  BE  PLOTTED  < 1 -YES,0-N0>?" ,Go_on 
IF  Go_on- 1  THEN  Repeat 
PRINT  "PU  SM" 

BEEP 

INPUT  "DO  YOU  LIKE  TO  DRAW  A  STRAIGHT  LINE?".Go_on 

IF  Go_on- 1  THEN 

BEEP 

INPUT  "ENTER  THE  SLOPE". SI 
BEEP 

INPUT  "ENTER  THE  INTERCEPT" , Ac 
FOR  Xa-Xmin  TO  Xmax  STEP  (Xmax-Xmin) 


Ya-Ac+Sl»Xa 
Y-(Ya-Ymin)«Sfy 
X-<Xa-Xffun>*$fx 
IF  Y<0  THEN 
Xan-< Ymin-Ac)/Sl 
X-TXa/n-Xmln) *Sf  x 
Y-0 


END  IF 

IF  Y> 1 00  THEN 
Xam*< Ymax-Ac)/SI 
X-<Xajn-Xmin)*Sfx 
Y-100 
END  IF 

PRINT  "PA" , X , Y , "PD" 
NEXT  Xa 
END  IF 

PRINT  "PU  SPO" 

END 
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10  •  FILE  NAME ;  Q_L0SS 

20  !  LAST  UPDATED:  23  MAY  33 

100  PRINTER  IS  701 

110  CLEAR  709 

120  BEEP 

130  INPUT  "ENTER  MONTH.  DATE  AND  TIME  <MM:DD:HH:MM:SS>" .DateS 
140  OUTPUT  709 : ”TD" :DateS 

150  OUTPUT  709 : "TD” 

180  ENTER  709;DateS 

170  PRINT  USING  "1  OX, ""Month,  date  and  tine:  1 4A" ;Date$ 

180  BEEP 

190  INPUT  "ENTER  INPUT  MODE  < 1 -3054A. 2-FILE)" , Im 
200  IF  Im-1  THEN 
210  BEEP 

220  INPUT  "GIVE  A  NAME  FOR  DATA  FIL£",D_f ileS 
230  CREATE  BDAT  D_file$,15 
240  ELSE 
250  BEEP 

260  INPUT  "GIVE  THE  NAME  OF  THE  EXISTING  DATA  FILE" ,D_f i Ie$ 

270  BEEP 

280  INPUT  "ENTER  THE  NUMBER  OF  RUNS  STORED". Nrun 
290  END  IF 

300  ASSIGN  SFile  TO  D_fileS 
310  BEEP 

320  INPUT  "GIVE  A  NAME  FOR  THE  OUTPUT  FILE" ,0_f i leS 

330  CREATE  BDAT  0_fileS.15 

340  ASSIGN  AFileo  TO  0_fileS 

350  K-0 

360  BEEP 

370  INPUT  "OK  TO  GO  AHEAD  (HIT  ENTER)". Ok 

380  IF  Im-1  THEN 

390  OUTPUT  709:"AR  AF19  AL21" 

400  OUTPUT  722: "FI  R1  T1  Z1  FL1” 

MO  FOR  1-0  TO  2 

420  OUTPUT  709: "AS  SA" 

430  ENTER  722:Eraf<I> 

440  IF  I>0  THEN  Emf (I)-ABS(Emf < I) >*10~6 

450  NEXT  I 

460  Ts  team»FNTvsv( (Emf(1)+Emf(2))*.5> 

470  OUTPUT  709 : "AR  AF23  AL24" 

480  OUTPUT  722; "FI  R1  T1  Zl  FL1  ** 

43U  h Ok  j-0  TO  1 

500  OUTPUT  709; "AS  SA" 

510  ENTER  722 ;Emf ( J+3) 

520  Enf  ( J+3) -ABS(Emf  <  J>3  > ) » 1 0  ~6 
530  NEXT  J 

540  GUTPUT  9Fiie;Emf<*) 

550  ELSE 

560  ENTER  5F  i  le:Emf  («) 

570  END  IF 

580  Ts team-FNTvsv( (Emf ( 1 ) ♦Emf (2) )* . 5) 

590  T  room*FNTvsv( (Emf ( 3) +Emf <  4 ) )» . 5) 

600  Power -(Emf  <0>*100)'2/5.76*1 .E-3 
610  Dt-Tsteam-Troom 

620  PRINT  USING  ”1  OX . 3D . 3D ,3(D0D . DD .2X>” : Power . Tsteam. Troom ,Dt 

630  K-K*1 

640  IF  Im-1  THEN 

650  BEEP 
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660 

670 

680 

690 

700 

710 

720 

730 

740 

750 


INPUT  "ARE  YOU  TAKING  MORE  DATA  ( 1-YES, 0-N0>?‘*.Go  .on 

IF  Go_on-1  THEN  380 

ELSE 

IF  K<Nrun  THEN  380 

END  IF 

END 

DEF  FNTvsv(V) 

T-V»< .0261 733441 6-V*<9 .2447859E-7-V-6 .  0746642E- 11)) 

RETURN  T 

FHEND 


Month,  date  and  tine:  06 : 0 1 :  1 4 : 55 :20 

The  following  analysis  was  performed  for  data  m  file  DD10 

Position  number  :  1  23456 

Wall  temperature  (Deg  C>  :  66.92  65.43  62.00  57.65  65.81  58.03 
Average  wall  temperature  ■  61.56  <Deg  C> 

Inlet  water  temperature  ■  23.94  (Deg  C> 

Outlet  water  temperature  •  26.32  (Deg  C> 

Log-mean- temp  difference  »  36.42  (Deg  C) 

Water  velocity  *  1.37  (m/S> 

Sieder-Tate  parameter  •  5342.8 

Nusselt  number  *  158.72 

Position  number  :  1  23456 

Wall  temperature  <Deg  C)  :  68.20  66.10  61.57  55.55  65.90  58.93 
Average  wall  temperature  *  61.42  (Deg  C) 

Inlet  water  temperature  ■  25.36  < Deg  C> 

Outlet  water  temperature  *  27.39  (Deg  C) 

Log-mean- temp  difference  ■  35.03  (Deg  C> 

Water  velocity  ■  1.92  (m/S) 

Sieder-Tate  parameter  ■  7032.7 

Nusselt  number  ■  199.76 

Position  number  :  1  23456 

Wall  temperature  (Deg  C)  :  68.10  65.85  60.61  54.03  65.49  57.68 
Average  wall  temperature  *  60.48  (Deg  C> 

Inlet  water  temperature  *  26.77  (Deg  C) 

Outlet  water  temperature  *  28.59  (Deg  C> 

Log-mean- temp  difference  »  32.79  (Deg  C> 

Water  velocity  *  2.47  (m/S) 

Sieder-Tate  parameter  *  8634.8 

Nusselt  number  *  247.76 

Position  number  :  1  23456 

Wall  temperature  (Deg  C)  :  67.70  65.11  59.70  52.58  64.40  56.49 
Average  wall  temperature  “  59.34  (Deg  C> 

Inlet  water  temperature.  ■  28.23  (Deg  C> 

Outlet  water  temperature  *  29.94  (Deg  C> 

Log-mean- temp  difference  *  30.22  (Deg  C) 

Water  velocity  ■  3.01  (m/S) 

Sieder-Tate  parameter  •  10181.6 

Nusselt  number  *  302.10 

Position  number  :  1  23456 

Wall  temperature  (Deg  C>  :  66.69  64.08  58.48  51.22  63.06  55.10 
Average  wall  temperature  *  58.04  (Deg  C) 

Inlet  water  temperature  ■  29.73  (Deg  C) 

Outlet  water  temperature  *  31.19  (Deg  C) 

Log-mean-temp  difference  *  27.58  (Deg  C) 

Water  velocity  *  3.56  (m/S) 

Sieder-Tate  parameter  *  11676.1 

Nusselt  number  ■  347.19 

Position  number  :  1  23456 

Wall  temperature  (Deg  C)  :  67.18  64.46  58.95  51.63  63.76  55.68 
Average  wall  temperature  *  58.53  (Deg  C) 

Inlet  water  temperature  *  31.91  (Deg  C) 

Outlet  water  temperature  *  33.28  (Deg  C) 
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Log-mean- temp  difference  *  25.93  (Deg  C) 

Water  velocity  *  3.84  (m/S> 

Sieder-Tate  parameter  ■  12518.8 

Nusselt  number  ■  371.81 

Position  number  :  1  2  3  4  5  6 

Wall  temperature  (Deg  C>  :  70.32  87.38  62.08  55.23  66.98  59.06 

Average  wall  temperature  *  61.84  (Deg  C> 

Inlet  water  temperature  *  33.57  (Deg  C) 

Outlet  uater  temperature  »  35.09  (Deg  C) 

Log-mean-temp  difference  *  27.51  (Deg  C> 

Water  velocity  »  3.29  (m/S> 

Sieder-Tate  parameter  »  11223.3 

Nusselt  number  *  331.04 

Position  number  :  1  23456 

Wall  temperature  (Deg  C>  :  74.00  71.28  65.83  58.91  70.63  62.82 

Average  ual 1  temperature  ■  65.70  (Deg  C> 

Inlet  water  temperature  *  35.46  (Deg  C> 

Outlet  uater  temperature  ■  37.16  (Deg  C) 

Log-mean- tenp  difference  *  29.38  (Deg  C> 

Water  velocity  ■  2.74  (m/S) 

Sieder-Tate  parameter  *  9850.0 

Nusselt  number  *  284.11 

Position  number  :  1  23456 

Wall  temperature  (Deg  C)  :  77.43  74.65  69.49  62.86  74.08  66.46 
Average  wall  temperature  *  69.42  (Deg  C) 

Inlet  uater  temperature  ■  37.27  (Deg  C> 

Outlet  uater  temperature  “  39.19  (Deg  C) 

Log-mean-temp  difference  -  31.18  (Deg  C> 

Water  velocity  *  2.19  (m/S) 

Sieder-Tate  parameter  *  8354.8 

Nusselt  number  -  240.23 

Position  number  :  1  23456 

Wall  temperature  (Deg  C)  :  80.70  77.96  73.09  67.18  77.57  70.37 
Average  uall  temperature  *  73.22  (Deg  C) 

Inlet  uater  temperature  *  38.92  (Deg  C) 

Outlet  uater  temperature  •  41.23  (Deg  C> 

Log-mean- temp  difference  *  33.12  (Deg  C> 

Water  velocity  *  1.64  (m/$> 

Sieder-Tate  parameter  -  6722.9 

Nusselt  number  ■  200.30 

Sieder-Tate  Coefficient  ■  .0294 

Least-Squares  Line: 

Slope  -  2.93533E-02 

Intercept  ■  0 . 00000E+00 
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Month,  date  and  time  : 05 : 30 : 1 5 : 27  : 1 5 

The  following  analysis  was  performed  for  data  i 

Measured  pressure  * 

Pressure  corresponding  to  measured  steam  temp  • 
Moncondensable  gas  concentration 
Noncondensable  gas  concentration  * 


Mater  velocity 
Reynolds  number 
Inside  heat-tran  coef 
Overal 1-heat-trans  coef 
Heat  flux 

Outside  heat-tran  coef 
Nusselt  coefficient 
Nus  coef  (from  heat  flux> 

Hater  velocity 
Reynolds  number 
Inside  heat-tran  coef 
Overal 1-heat-trans  coef 
Heat  flux 

Outside  heat-tran  coef 
Nusselt  coefficient 
Nus  coef  (from  heat,  flux) 

Hater  velocity 
Reynolds  number 
Inside  heat-tran  coef 
OveTal 1-heat-trans  coef 
Heat  flux 

Outside  heat-tran  coef 
Nusselt  coefficient 
Nus  coef  (from  heat  flux) 

Hater  velocity 
Reynolds  number 
lnsioe  heat-tran  coef 
Overal l-heat-trans  coef 
Heat  flux 

Outside  heat-tran  coef 
Nusselt  coefficient 
Nus  coef  (from  heat  flux) 

Hater  velocity 
Reynolds  number 
Inside  heat-tran  coef 
Overal 1-heat-trans  coef 
Heat  flux 

Outside  heat-tran  coef 
Nusselt  coefficient 
Nus  coef  (from  heat  flux) 


-  0.82  (m/S) 

-  13139.2 

-  5356.2  (W/m  2-0 

■  3393.6  <W/m'2-K> 

-  2 . 22899E+05  <H/m'2> 

-  19990.5  (W/m"2-K) 

-  1 1 304 .5  (H/m'2-K) 

-  10797.9  (W/m'2-K) 

-  1 .64  (m/S) 

-  26418.9 

-  9074.3  (H/m'2-K) 

-  4109.5  ( W/m"2-K) 

-  2.68764E+05  (H/m'2) 

-  11186.9  (H/m~2-K> 

-  10320.1  (H/m'2-K) 

-  10005.0  (W/nT2-IO 

-  2.47  (m/S) 

-  40022.4 

-  12466.1  (H/m'2-K) 

-  47?5.9  (H/n'2-K) 

-  3.06S51E+05  (H/m'2) 

-  10550.9  < W/m'2-K) 

-  9822.2  (H/m'Z-K) 

-  9484.8  (W/m^-K) 

-  3.29  'm/S> 

-  53980.1 

■  15667.3  (H/m‘2-K) 

-  5331 .5  <W/m'2-K> 

-  3 . 43200E+05  <U/nT2> 

-  10926.6  (W/ra~2-K> 

-  95M.1  <M/m‘2-K) 

-  9075.4  (W/n'2-K) 

-  3.84  <m/S> 

■  63586.5 

-  17728.9  <W/m'2-K) 

-  5652.2  (W/m‘2-K) 

-  3.61307E+05  (H/m'2) 

-  11060.4  <W/m"2-K) 

-  9359.5  (H/n'2-K) 

-  8889.1  <H/m"2-K) 


NOTE:  5  X-Y  pairs  were  stored  in  plot  data  fil 


file  DR1 

9 . 578E+04  (Pa) 
9.010E+04  (Pa) 
5.8  (No  1  /.) 
9.0  ( Mass  V. ) 


PDR1 
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END  IF 

20d0 

PRINT 
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2090 

PRINT 

MMt 

USING 

"12X, 

2100 

PRINT 

USING 
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2110 

PRINT 

USING 

"12X, 
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'“‘Wall  temperature  (Deg  C) 
''"Average  wail  temperature 
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2 


4 


'"'  ,6<D0  .DD. IX)": Tu<«> 
'“'.DO. DO.""  (Deg  C>"" 


;  Tual 1 


2650  FNEND 

2H0  Mu-K57609473E-3-T*<3.51  1 98576E-5-T«<  3. 583581 6E-?-T»1 .3655861  t5E-9>> 

2680  RETURN  Hu 

2690  FNEND 

2700  DEF  FNKy<T) 


Position  number  :  1  2  3  A  5  6 

Mall  temperature  (Deg  C)  :  67. !8  64. A6  58.95  51.63  63.76  55.68 
Average  wall  temperature  *  58.53  (Deg  C) 

Inlet  water  temperature  *  3t.St  (Deg  C> 

Outlet  water  temperature  ■  33.28  (Deg  C) 
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